PROCEEDINGS 

OF SCIENCE 



Event shapes at hadron colliders 



Andrea Banfi* 1 

ETH Zurich, Institute for Theoretical Physics 



E-mail: banfi@itp. phys . ethz . ch 



We present precision results for distributions in global event shapes that can be measured at 
hadron colliders within experimental limitations. These predictions are obtained by combining 
exact next-to-leading order (NLO) with the all-order resummation of large logarithms of soft and 
collinear origin. We then discuss how event-shape measurements can be used for the tuning of 
Monte Carlo event generators, for tests of models of hadronisation and underlying event, and as 
discriminatory tools between QCD jet-like and New Physics events. 
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1. Introduction 

Event-shape variables are infrared and collinear safe measures of the geometrical properties 
of the hadronic energy momentum flow, giving an idea on whether an event is pencil-like, planar, 
spherical, etc. Measurements of their mean values and distributions have played a crucial role at 
LEP, for precise determinations of the strong coupling a s , for tests of analytical models of hadro- 
nisation corrections, and for validation of Monte Carlo (MC) event generators (see ^ for a recent 
review). In spite of the success of these studies in e + e~ annihilation, very little attention has been 
devoted to their counterparts at hadron colliders [^, [|]. This was mainly because the quantities 
that were conveniently measured experimentally could not be accurately computed in perturbative 
QCD. Here I would like to present event shapes that can be measured at actual hadron colliders, 
and whose distributions can be computed in perturbative QCD at the accuracy needed to have a 
reliable estimate of the associated theoretical uncertainties. 

2. Event-shape variables at hadron colliders 

In hadron-hadron collisions, we consider events with two hard central jets, and define event 
shapes that vanish in the limit of two narrow jets. For instance, given all hadrons {qf\ in a rapidity 
region ^ (for instance |tj,-| < r]^) and using transverse momenta q±j only, we define the transverse 
thrust 7j_ and the thrust minor T m as follows 

r ± =max%%^, r w = £ ' eyl *t' X ? xl - (2.D 

«x Li^\<lX,i\ Lie^\l-L,i\ 

We can also introduce boost-invariant event shapes involving longitudinal degrees of freedom, like 
invariant masses or broadenings, or the three-jet resolution parameter y$. An extensive list of 
hadronic event-shape definitions can be found in [Q, 

There are three basic reasons why dijet event shapes can be studied experimentally with very 
first data. First, cross sections for dijet production are large both at the Tevatron and at the LHC, 
as shown in table [j] From table [I] one can also see that the flavour content of a sample can be 
varied by changing the leading jet jo r cut. Low-p f samples (Tevatron with p t \ > 50GeV, LHC 
with pti > 200GeV) are gluon dominated, while initial-state quarks become more important for 
high-p t samples (Tevatron with p t \ > 200GeV, LHC with p t \ > ITeV). Second, event shapes are 





LO NLO 


qq->qq qg^qg gg^gg 


Tevatton, p t \ > 50GeV 
Tevatron, p n > 200GeV 
14TeV LHC, p n > 200GeV 
14TeV LHC, p n > ITeV 


60nb 116nb 
59pb lOlpb 
13.3nb 23.8nb 
6.4pb 10.5pb 


10% 43% 45% 
41% 43% 12% 
7% 40% 50% 
31% 51% 17% 



Table 1: Cross sections for the production of two jets in a central rapidity region ( |_yj ets | < 0.7 at the Tevatron 
and \y-ps | < 1 at the LHC) with a cut on p t \, the transverse momentum of the leading jet. On the right it is 
possible to see the relative importance of each partonic subprocess. 

normalised quantities: experimental uncertainties associated with jet energy scale cancel between 
numerators and denominators, see eq. (|Q|). Finally, since one usually measures normalised differ- 
ential distributions, like 1 / o do / dT m , no determination of luminosity is required. 
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From a theoretical point of view, event-shape distributions can be computed at NLO with NLO- 
JET++ However, as can be seen in the example plots of figure [j], both LO and NLO predictions 
diverge at small values of T„ jg (the directly global variant of T m , see section ||). Only a combi- 
nation of NLO and next-to-leading logarithmic (NLL) resummation (labelled NLL+NLO) gives 
a distribution that is sensible for any value of T m g . In particular, resummation restores the cor- 
rect physical behaviour at T m>g — >• 0, corresponding to vanishing probability of having accelerated 

NLL resummation involves writing the 
integrated T m>g distribution as an exponent 
Z(T ltl:g ) = exv[Lgi(a s L) + g 2 (a s L)], with 
L = ln(l/r m ? ), gi (a s L) resumming the lead- 
ing logarithms (LL, a"L n+l ), and g 2 (a s L) 
the NLL, a"L". Knowledge of g\ (a s L) de- 
termines the position of the peak of the dif- 
ferential distribution, typically in the region 
a s L ~ 1 . In the peak region g 2 (cc s L) becomes 
of order one, and is therefore needed to sta- 
bilise both the position and the height of the 
peak. 

NLL resummability is guaranteed for 
variables satisfying the following condi- 
tions [|7p: a specific functional dependence 
on a single soft and collinear emission; (con- 
tinuous) globalness, i.e. sensitivity to emis- 
sions in the whole of the phase space; recursive infrared and collinear (rlRC) safety, a subtle 
mathematical condition on the event-shape scaling properties with multiple emissions. If these 
conditions are satisfied, the relevant emissions that contribute to event-shape distributions at NLL 
accuracy are soft and collinear parton clusters widely separated in rapidity, which, due to QCD 
coherence, can be considered as independently emitted from hard legs. Since a similar pattern of 
emissions is simulated by MC event generators, one may expect that most features of rlRC safe 
global event-shape distributions are correctly described by these theoretical tools. 

Of the three constraints, the most difficult to satisfy experimentally is globalness, due to the 
fact that the measurement region is preferably restricted to the central detector region (e.g. \r[ | < 
2.5 at the LHC), and in any case no measurement is actually performed in the very forward regions 
(corresponding to a limiting rapidity r\ c = 5 at the LHC). However [Q], one can devise classes of 
global event shapes even at hadron colliders: directly global, where the region ^£ extends up to 
the maximum available rapidity t] c ; exponentially suppressed, where ^ is inside the acceptance 
of the central detectors (e.g. |t]| < 1 at the Tevatron and |t]| < 1.5 at the LHC), while outside this 
region we add to the event-shape definition a term that exponentially suppresses the contribution 
of hadrons in the forward regions; recoil, where measurements are performed only in a central 
region % ', and we add a term that is sensitive to emissions outside ^ through recoil. In the last case 
however a numerical breakdown of NLL resummation occurs in the region where the event shape 
is small. 
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3. Event-shape distributions the Tevatron and at the LHC 

In []|] we have performed a NLL+NLO 
resummation for a number of selected event 
shapes. Figure |2| shows predictions for the di- 
rectly global thrust minor T nug , together with 
theoretical uncertainties. The latter, aimed at 
having an indication of missing NNLO and 
NNLL corrections, are estimated via: asym- 
metric variation of renormalisation and factori- 
sation scales jXR and p F ; variation of the log- 
arithm to be resummed \n(XT mg ) with 1/2 < 
X < 2; variation of the matching scheme (log- 
R or mod-R). We observe that uncertainties are 
under control and within ±20% in a wide range 
of values of T„ hg . Similar results are obtained 
for all considered event shapes. 

These predictions are valid at parton level 
only, so it is interesting to investigate the impact 
of hadronisation and underlying event on event- 
shape distributions. This can be done with 
MC event generators. Figure || shows that k t - 
algorithm jet resolution parameters, for exam- 
ple j3 g , are essentially not affected by hadroni- 
sation and underlying event, while event shapes, 

like Pt.s, get moderate hadronisation corrections, falling as an inverse power of the jet p t , but get 
a huge contribution from the underlying event. This different sensitivity shows that event-shape 
distributions can be exploited for the validation of MC event generators. Jet resolution parameters 
are better suited for tunings of parton shower parameters, while with remaining event shapes one 
can test models of hadronisation and underlying event. 




Figure 2: NLL+NLO predictions and theoretical 



uncertainties for the distribution in T, 
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Figure 3: Distributions in directly global three-jet resolution „ (left) and exponentially suppressed total 
mass p T g (right), as obtained with the MC event generator PYTHIA [||] for LHC with = 14TeV. 
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4. Event shapes for New Physics 

A common use of event shapes is that of discriminating among events with different topolo- 
gies. This is particularly important in New Physics searches, where one expects events with mas- 
sive particles to be much broader than dijet events. We have then tried to assess the performance 
of known hadronic event shapes for such studies. First, when considering symmetric events with 
an arbitrary number of particles in the transverse plane, one can only distinguish between two- and 
multi-jet events, irrespectively of the number of jets. One can then try to discriminate among dif- 
ferent topologies in a sample with the same number of jets (three jets in the considered case). One 
then finds (see figure |]) that infrared and collinear safe variables fare much better in this respect 
than unsafe ones (like the widely used transverse sphericity). We remark also that event shapes like 




0.2 0.4 0.6 0.8 1 1.2 0.2 0.4 0.6 0.8 1 1 .2 0.2 0.4 0.6 0.8 1 1.2 



Figure 4: Distribution in the value of three different variables obtained from events with two different hard 
configurations dressed with parton shower. The two starting configurations are a generic three-parton and a 
totally symmetric (Mercedes) event in the transverse plane. 

the broadenings, which treat tranverse and longitudinal degrees of freedom on equal footing, are 
better suited for the identification of massive particle decays, since their value is hardly affected by 
the orientation of the event plane. 

For practical applications, it is however desirable to have variables that are more sensitive to the 
spherical limit. One example presented in ^ is the supersphero observable which is non-zero only 
for events in which there are three non-coplanar particles in each of the "hemispheres" in which 
the event is divided by the transverse thrust axis. We believe that phenomenological applications 
of variables like supersphero, as well as better final-state observables for New Physics, constitute 
an important subject that deserves further studies. 
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